Coherent THz emission of intrinsic Josephson junction stacks in the hot spot regime 
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We report on THz emission measurements and low temperature scanning laser imaging of 
Bi2Sr2CaCu2 0g intrinsic Josephson junction stacks. Coherent emission is observed at large dc 
input power, where a hot spot and a standing wave, formed in the "cold" part of the stack, coexist. 
By varying the hot spot size the cavity resonance frequency and the emitted radiation can be tuned. 
The linewidth of radiation is much smaller than expected from the quality factor of the cavity mode 
excited. Thus, an additional mechanism of synchronization seems to play a role, possibly arising 
from nonequilibrium processes at the hot spot edge. 

PACS numbers: 74.50.+r, 74.72.-h, 85.25.Cp 



Phase synchronization is one of the prerequisites to 
use Josephson junction arrays as tunable high frequency 
sources [l| . While Nb based junctions are limited to fre- 
quencies well below 1 THz, intrinsic Josephson junctions 
(IJJs)[i,i] in BiaSraCaCuaOg (BSCCO) are, at least in 
principle, able to operate up to several THz. Stacks of 
many junctions can be made e.g. by patterning mesa 
structures on top of single crystals. For many years, 
investigations focused on small structures consisting of 
some 10 IJJs, with lateral mesa sizes of a few /im. Here, 
with few exceptions [3, la] , the IJJs in the stack tended 
to oscillate out-of-phase or were not synchronized at all. 
Experimental and theoretical studies included the gen- 
eration of synchronous Josephson oscillations by moving 
Josephson vortices |4l-[l0|. the use of shunting elements 
[ll| - ll4| . the excitation of Josephson plasma oscillations 
via heavy quasiparticle injection |15l . |l6[ or the investiga- 
tion of stimulated emission due to quantum cascade pro- 
cesses jl7| . High-frequency emission of unsynchronized 
intrinsic junctions has been observed up to 0.5 THz jlSj . 

Recently, coherent off-chip THz radiation with an ex- 
trapolated output power of some /^W was observed from 
stacks of more than 600 IJJs, with lateral dimensions 
in the 100 /xm range |19l - [23| . Phase synchronization in- 
volved a cavity resonance oscillating along the short side 
of the mesa. This radiation was studied theoretically in 
a series of recent papers, either based on vortex- type or 
plasmonic excitations [2J-|33| , or on nonequilibrium ef- 
fects caused by quasiparticle injection |34| . 

While in experiments |19l - l23| THz emission was ob- 
tained at relatively low bias currents and moderate dc 
power input, using low temperature scanning laser mi- 
croscopy (LTSLM) we have shown that standing wave 
patterns, presumably associated with THz radiation, can 
be obtained at high input power where, in addition, a hot 
spot (i.e. a region heated to above the critical temper- 
ature Tc) forms within the mesa structure j35|. Hot spot 
and waves seem to be correlated. The purpose of the 
present work is to investigate THz emission in this high 



power regime in detail, combining THz emission measure- 
ments and LTSLM. Apart from further clarifying the role 
of the hotspot, we are specifically interested in the ques- 
tion whether coherent radiation can be achieved in spite 
of the high temperatures involved. 

For the experiments BSCCO single crystals were grown 
by the floating zone technique in a four lamp arc-imaging 
furnace. Below we discuss results from two samples. 
Sample 1 was patterned on a crystal that was annealed 
in vacuum at 650 °C for 65 hours. It had a Tc of 86.6 K 
and a transition width ATc of 1.5 K. Sample 2 was made 
on a crystal {Tc = 87.6 K, ATc = 1.5 K) annealed in vac- 
uum at 600 °C for 72 hours. To provide good electrical 
contact the single crystals were cleaved in vacuum and 
a 30 nm Au layer was evaporated. Then, conventional 
photolithography was used to define the mesa size in the 
a-b plane (length 330 /xm, width 50 /xm for both samples). 
Ar ion milling yielded mesas with measured thicknesses 
of, respectively 1 ^m (sample 1) and 0.7 /im (sample 2) 
along the c-axis (corresponding to, respectively, stacks of 
670 and 470 IJJs). Insulating polyimide was used to sur- 
round the mesa edge at which a Au wire was attached to 
the mesa by silver paste. Other Au wires were connected 
to the big single crystal pedestal as grounds. In order 
to provide a load line for stable operation, the mesas 
were biased using a current source and variable resistor 
in parallel to the mesa, cf. Fig. 1 in Ref. |35|. The volt- 



age measured across the mesa includes the resistance of 
the contacting Au wires and the resistance between these 
wires and the mesa (4ri for sample 1, 6.5ri for sample 2). 
In the data discussed below this resistance is subtracted. 
THz emission measurements were performed in 
Tsukuba; the samples were subsequently shipped to 
Tiibingen for LTSLM. In total we detected THz emis- 
sion from 7 (out of a total of 12) mesas on 5 different 
crystals. The LTSLM setup is described in Ref. j3a| . In 
brief, the beam of a diode laser (modulated at 10-80 kHz, 
spot size 1-2 jimj is deflected by a scanning unit and fo- 
cused onto the sample surface. Local heating by 2-3 K in 
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Figure 1: (color online). Schematic of the THz emission setup 
together with a sketch of the sample. 



an area of a few /zm^ and about 0.5 iim. in depth causes 
a change AV (detected using lock-in techniques) of the 
voltage V across the mesa serving as the contrast for the 
LTSLM image. Standing waves can be imaged due to the 
beam-induced local change of the quality factor, leading 
to a strong signal AV at antinodcs and a weak signal at 
nodes. The edge of a hot spot leads to a strong signal, 
while AV is small in the interior of the hot spot and in 
the " cold" part of the mesa [3^ . The THz emission setup 
is shown schematically in Fig[T] together with a sketch of 
the s amp le. The interferometer is similar to the one in 
Ref. l36l The crystal is mounted in a continuous flow 
cryostat with a polyethylene window. The emitted radi- 
ation, chopped with a frequency of 13 Hz, is collected by 
a 90° off-axis parabolic mirror (2" diameter, focal length 
100 mm) and deflected towards two lamellar split mirrors 
dividing the incoming radiation into two beams with al- 
most the same intensity. A phase difference is introduced 
by moving one of these mirrors with a translation stage. 
Via a second parabolic mirror the radiation is deflected 
to the bolometer (gauged by using a Gunn oscillator for 
reference) . A 1 THz cut-off type filter is used in front of 
the bolometer. The numbers we quote below are based 
on the detected power, ignoring additional losses (by a 
factor of 2-5) from the continuous flow cryostat window 
and the mirrors. The solid angle of our setup, defined by 
the aperture of the Winston cone in front of the bolome- 
ter, is 0.04 sr. The frequency resolution, depending on 
the maximum mirror displacement, is 10-15 GHz. 

We first discuss THz emission data of sample 1, mea- 
sured at T =20 K. Figure [5] (a) shows by open circles 
the current voltage characteristic (IVC), as measured in 
Tsukuba. Starting from zero current, groups of junc- 
tions become resistive for / > 20mA. At / = 30 mA 
the whole stack switches to the resistive state. Lower- 
ing the bias in this state the IVC is continuous down to 
about 10 mA where another switch to a current of about 
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Figure 2: (color online). IVCs (a) of sample 1 at T = 20K, 
as measured in Tsukuba [open (red) circles] and in Tiibin- 
gen [solid (blue) circles]. Labelled arrows refer to bias values 
where LTSLM images are taken (c.f. Fig. |3]). Emission power 
detected by the bolometer as a function of current (b) and 
voltage (c). Black lines refer to high bias (/ >10 mA), grey 
(orange) lines to low bias (/ <7 mA). Arrows in (c) indicate 
where Fourier spectra (d) have been taken; the peak at 0.669 
THz corresponds to 1/ = 944 mV and I =17.2 mA, the peak 
at 0.753 THz to 1/ = 1046 mV and / =13.2 mA. 



8 mA occurs. As revealed by LTSLM this switch corre- 
sponds to the disappearance of the hot spot present at 
high bias. For 3 mA < / < 10 mA the IVC is continuous 
until further switches and hystereses appear due to the 
retrapping to the zero voltage state of some of the I J Js in 
the stack. Figures [2] (b) and (c) show the emission power 
detected by the bolometer as a function of, respectively, 
/ and V. The largest response occurs in the high power 
(hot spot) regime and covers a current range from 10 mA 
up to about 30 mA. Here, the voltage across the mesa 
ranges from 0.83 to 1.13 V. If all N IJJs in the stack are 
frequency locked and participate in radiation we expect 
that the Josephson relation V/N = $o./ holds, where 
$0 is the flux quantum and / is the emission frequency. 
The Fourier spectrum, at given bias, revealed a single 
sharp emission line (with a resolution limited A/ « 12 
GHz lincwidth) continuously shifting with V. Two spec- 
tra are shown in Fig. [2] (d). From these data we infer 
N = 676 ± 5 in very good agreement with the 1 fj.m thick- 
ness of the mesa. The maximum detected power was 6 
nW, corresponding to roughly 2 ^W when extrapolating 
to 47r, even without taking additional losses into account. 
Thus the emitted power is comparable to the one quoted 



(a) 



■AV(nV) 
100 




(b) 





-inTro ' ' ' 


L 


75 


1 /I — ® 


- 


50 


- // 1 


- 


25 


-ff \ 


- 


n 


/ 


\ 



300 x(nm) 



100 200 
x(ijm) 



300 



Figure 3: (color online) LTSLM imabges (a) and linescans 
along the center of the mesa (b) of sample 1 at T = 20 K, 
corresponding to the bias points (i),(ii) a-nd (iii) indicated in 
Fig. [21^ a). In (a), image (ii b) shows a simulated standing 
wave pattern superposed on the data of image (ii). 



in Ref. [lOl Sample 1 also showed a bolometric response 
at low bias [c.f. grey (orange) lines in Figs. [H (b) and 
(c)]. Here, the detected power (of order 0.4 nW) was too 
low for Fourier spectrometry. 

The different peaks visible in the radiometer data of 
Figs. [2] (b) and (c) are likely to correspond to different 
cavity modes excited in the mesa (a similar case is shown 
in Fig. 3 of Ref. [Sa]). For the present sample, only for a 
bias near the main emission peak the standing wave pat- 
tern appeared clear enough to be analyzable. Fig. [3l^a) 
shows LTSLM images at the bias points (i), (ii) and (iii) 
indicated on the IVC shown in Fig. [5] (a) by solid cir- 
cles. Fig. [3] (b) shows the corresponding linescans, taken 
along the center of the mesa. The bright feature visi- 
ble in the left part of the mesa is the hot spot growing 
with increasing input power. Image (ii) shows a stand- 
ing wave pattern located in the right part of the mesa. 
There are 3 half waves along the long side of the mesa 
and one half wave along the short side. A simulation of 
this wave pattern, superposed to the data at bias point 
(ii), is shown in image (ii b). Here we assumed that the 
LTSLM signal is dominantly sensitive to the magnetic 
field B, i. e. A^ oc B^. From the LTSLM data we infer 
Aj: « Ay « 100 fim. Using / — c(A~^ -I- A~^ 



)i/2, with / 
= 0.67 THz, we find a mode velocity c w 4.7 • lO^m/s, 
which is fully compatible with the in-phase mode veloc- 
ity at the elevated temperatures (60-70 K) estimated for 
the cold" part of the mesa [3^, |33]- Other (collective) 
mode velocities are by a factor of at least 2 lower. If the 
junctions would oscillate incoherently, the relevant mode 
velocity is the Swihart velocity c « 3 • lO^m/s, which is 
even lower. Wc thus conclude that all junctions oscillate 
in-phase. 

We next discuss emission data of sample 2. IVCs, the 
bolometric response and Fourier spectrometer data are 
shown in Fig. [4] The largest emission was detected near 
a voltage of 456 mV, c.f. Fig. 1 (b), (c). Fig. I^d) 
shows Fourier spectra at V =440 mV and V =451 mV. 
From the emission peaks at, respectively, 531 GHz and 
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Figure 4: (color online). IVCs (a) sample 2 at T = 40 K. 
Emission power detected by the bolometer as a function of 
current (b) and voltage (c). (d) shows two spectra of the 
emitted power at V = 440 mV and V = 451 mV. 
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Figure 5: (color online) Range of emission frequencies / (a) 
and detected emission power (b) of sample 2 as a function of 
bath temperature. Inset in (a) shows the voltage across the 
mesa as a function of /. Solid line in the inset is function 
V ^ Nf^o, with iV=400. 



545 GHz wc find the Joscphson relation V/N = $0/ 
to hold, with N « 400 {i.e. somewhat less than the 
470 junctions expected from the meas thickness). The 
16 GHz linewidth is at the resolution limit. LTSLM re- 
vealed that a hot spot formed for / > 25mA, covering 
about half of the mesa area for the bias where maxi- 
mum emission was detected. A standing wave pattern 
was vaguely visible, however further analysis of the LT- 
SEM data was prevented by a spurious strong response. 



presumably caused by additional IJ Js switching to the re- 
sistive state when illuminated by the laser. Fig. [S] shows 
emission data taken at different bath temperatures. The 
frequency range of detectable emission is shown in Fig. 
[5] (a). Radiation was found between 15 K and 58 K. 
At the lower temperatures the emission frequencies var- 
ied continuously between 600 GHz and about 750 GHz; 
with increasing bath temperature the frequency of max- 
imum emission power decreased, to about 450 GHz near 
55 K; also, the frequency tunability decreased with in- 
creasing bath tempeature but still was on the order of 
10% at 50 K. In all cases we found the Josephson re- 
lation V/N = $o/ to hold [c.f. inset of Fig. [5] (a)], 
with N = 400. Thus, in contrast to the "low bias" case, 
where the emission frequency is essentially fixed by the 
geometry of the mesa, it can be tuned in the hot spot 
regime. This can be understood if one accepts that one 
of the roles of the hot spot is to vary the length of the 
"cold" part of the mesa until, eventually, resonant con- 
ditions are achieved. Further, both with increasing dc 
input power and with increasing bath temperature, the 
(in-phase) mode velocity decreases, from about 7-lO^m/s 
for low T, towards zero when the cold part of the mesa 
reaches T^ [35|- Thus, the change of / with bath temper- 
ature essentially reflects the decrease in mode velocity. 

Still, the hot spot might even take an active role in 
THz generation and/or synchronization. For sample 1, 
the (resolution limited) 12 GHz linewidth of radiation 
is more than 50 times smaller than the oscillation fre- 
quency. For sample 2, A/// < 0.03 at T = 40 K. If, in 
the hot spot regime. A/ were determined by the quality 
factor Q of the cavity, one would require a Q = f / ^f of 



at least 30-50, which is very unlikely to be the case at 
the elevated temperatures of the hot spot regime, where 
we expect Q < 5. Thus, although the cavity seems to be 
a necessary ingredient for synchronization, the (nonequi- 
librium) phenomena, in the spirit of Rcfs. J15l4l7l |3J], 
taking place at the hot spot edge may play an impor- 
tant role as well in synchronizing the junctions, eventu- 
ally leading to the small values of A/ observed. The 
fact that the Josephson frequency- voltage relation holds 
rules out a purely quasiparticle related process, however, 
should not contradict a scenario of interacting Josephson 
and quasiparticle currents. 

In summary, we have investigated THz emissison from 
intrinsic Josephson junction stacks by off-chip emission 
measurements and by low temperature scanning laser 
imaging. The stacks emit radiation coherently at high 
input power where a hot spot has formed. Here, standing 
wave patterns are observed pointing to the importance of 
cavity resonances for synchronization. By changing the 
size of the "cold" part of the mesa the hot spot effec- 
tively adjusts the size of the cavity, making the emission 
frequency tunable. The linewidth of radiation, however, 
is too small to be determined by the quality factor of the 
cavity alone. An additional mechanism of synchroniza- 
tion seems to play a role and may arise from nonequilib- 
rium processes at the edge of the hot spot. 
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